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Abstract 



There is experimental and theoretical evidence that the broad rapidity distribution 
of net proton yield {p — p) in central heavy-ion collisions at SPS energies could 
be a signal of non-equilibrium properties of the system. We show that the broad 
rapidity shape can be well reproduced in the framework of a non-linear relativistic 
Fokker-Planck dynamics which incorporates non-extensive statistics and anomalous 
■ diffusion. 
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^ ' 1 Introduction 

Qh! In relativistic heavy-ion collisions nuclear matter is believed to be compressed 

r-| ! to high baryon density and nuclei undergoing central collisions reduce strongly 

their original longitudinal momentum. This loss of rapidity, usually referred 



as baryon stopping, is an important characteristic to understand the reaction 
mechanism at high energy density. 

Isotropic thermal models fail to reproduce the broad and flat rapidity distribu- 
tion observed in the central Pb-Pb collisions at 158 A GeV/c [1]. Collective and 
anisotropic flows must be considered (with the introduction of new parame- 
ters) to take into account the experimental spectra [2,3]. Complicate evolution 
codes take only partially into account of such distributions and the physical 
evolution of the system is not completely understood [4]. 

In this paper we study the evolution of the rapidity distribution from a macro- 
scopic point of view by using a non-linear relativistic Fokker-Planck equation 
and we show that the observed broad rapidity shape could be a signal of non- 
equilibrium properties of the system. Let us note that a similar approach, 
within a linear Fokker-Planck equation, has been previously studied in Ref . [5] 
by using a linear drift in the space of the rapidity and a free parameter diffu- 
sion coefficient. With this choice, the author found a strongly violation of the 



Preprint submitted to Elsevier Preprint 



1 February 2008 



fluctTiation-dissipation theorem. We will see that generalizing the Brownian 
motion to the relativistic kinetic variables, the standard Einstein relation is 
satisfied and Tsallis non-extensive statistics emerges in a natural way from 
the non-linearity of the Fokker-Planck equation. 



2 Non-linecir Fokker-Planck equation in relativistic dynamics 

A class of anomalous diffusions are currently described through the non-linear 
Fokker-Planck equation (NLFFE) 
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where D and J are the diffusion and drift coefficients, respectively. Tsallis and 
Bukman [6] have shown that, for linear drift, the time dependence solution of 
the above equation is a Tsallis-like distribution with q = 1 + — z/. The norm 
of the distribution is conserved for all times only ii /i — 1, therefore we will 
limit the discussion to the case u — 2 — q. 

We are going to study the time evolution of the particle distribution by means 
the above NLFFE in the rapidity space. Basic assumption of our analysis is 
that the rapidity distribution y is not appreciably influenced by transverse 
dynamics which is considered in thermal equilibrium. Such hypothesis is well 
confirmed by the experimental data and adopted in many theoretical works 
[1-3]. A crucial role in the solution of the above NLFFE plays the choice of the 
diffusion and the drift coefficients. Such a choice influences the time evolution 
of the system and its equilibrium distribution. 

Since the temperature at freeze-out exceeds 100 MeV, Boltzmann approxi- 
mation is usually adopted and the single particle equilibrium distribution is 
written as 



(X E exp(— £'/T) = m± cosh(|/) exp(— m_L cosh.{y)/T) , (2) 



where y is the rapidity, = y m'^ + p]_ is the transverse mass, T is the tem- 
perature. Imposing the validity of the Einstein relation for Brownian particle, 
we can generalize to the relativistic case the standard expressions of diffusion 
and drift coefficients as follows 

D^-fT, J(|/) = 7m±sinh(|/) = 7P|| , (3) 



where py is the longitudinal momentum and r is a common constant. It is easy 
to see that the above coefficients give us the Boltzmann stationary distribution 
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(2) in the linear Fokker-Planck equation {q = = while the stationary 
solution of the NLFPE (1) with i/ = 2 — g is a Tsallis-like distribution with 
the relativistic energy E = m_|_ cosh(?/). In the light of the above results, we 
generalize the results of Tsallis-Bukman by searching as a solution of the 
NLFPE (1) the following time dependent Tsallis-like distribution 

fgiy,m^,t) = {l- (1 - q) Pit) m^cosh[y - y^{t)]y^'^^~'^ . (4) 



The unknown functions = 1/T{t) and ym{t) have been derived by means 
numerical integration of Eq.(l) with initial 5-function condition depending on 
the value of the experimental projectile rapidities. The rapidity distribution 
at fixed time is then obtained by numerical integration over the transverse 
mass m_L (or transverse momentum) as follows 

dN 



dy 



{y, t) = c cosh(?/) /^(y, m_L, t) dm^ , (5) 
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where c is the normalization constant fixed by the total number of the par- 
ticles. Rapidity spectra calculated from (5) will ultimately depend on two 
parameters: the "interaction" time t = '-ft and the non-extensive parameter 
q. Therefore, no more free parameters are used in this analysis respect to 
previous theoretical studies [2-4]. 

3 Results and Conclusions 

We have studied the time dependent rapidity distribution (5) in the case of 
central Pb+Pb collisions at 158 GeV/c and our results are compared with the 
net proton experimental NA49 data [1]. The obtained spectra are normalized 
to 164 protons and the beam rapidity is fixed to j/cm = 2.9 (in the cm. 
frame) [1]. In Fig.l we show the calculated rapidity spectra compared with 
the experimental data. The full line corresponds to the NLFPE solution (5) 
at r = 0.82 and q = 1.25; the dashed line corresponds to the solution of the 
linear case (g = 1) at r = 1.2. Only in the non-linear case [q ^ 1) exists a 
(finite) time for which the obtained rapidity spectra well reproduces the broad 
experimental shape. A value of g 7^ 1 implies anomalous superdiffussion in the 
rapidity space, i.e., \y{t) — yuif)^ scale like with a > 1 [6]. 

We have studied the rapidity spectra from a macroscopic point of view. A 
good agreement with the experimental data is reached assuming a deviation 
from the equilibrium and a non-linear evolution which implies non-extensive 
Tsallis-like distribution. Non-extensive features in relativistic collisions are 

^ This result cannot be obtained if one assumes a linear drift coefficient as in Ref.[5]. 
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Fig. 1. Rapidity spectra for net proton production (p—p) in central Pb+Pb collisions 
at 158A GeV/c (grey circles are data reflected about ycm = 0) [!]• PuH line corre- 
sponds to our results by using a non-linear evolution equation (q = 1.25), dashed 
line corresponds to the linear case {q = 1). 

outlined in several works [7-11], even if a microscopic justification of these 
effects is still lacking. In Ref. [4] , the authors show that diquark breaking during 
inelastic collisions in dual parton model could be responsible for the observed 
baryon stopping. The possible connection between the nature of these inelastic 
collisions and the above observed non-extensive statistical effects is under 
investigations. 
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